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I summarize the derivation of a set of Feynman rules appropriate for the perturbative descrip- 
tion of the nonequilibrium dynamics of the symmetry-breaking phase transition in Xcj> 4 theory in a 
Robertson- Walker universe. The approximation scheme I develop provides for a treatment of dis- 
sipative effects which are essential to an adequate description of the nonequilibrium state. It also 
provides for the emergence from an initially symmetric state of a final state exhibiting the proper- 
ties of spontaneous symmetry breaking without the introduction by hand of any explicit symmetry 
breaking. 
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I. INTRODUCTION 

Thermal field theory strongly suggests that the hot, 
dense matter present in the rapidly expanding early uni- 
verse would undergo a variety of phase transitions. In 
particular, the "new" inflationary scenario of Linde || 
and Albrecht and Steinhardt envisaged a symmetry- 
breaking phase transition in a grand unified theory (mod- 
ifying the earlier, seminal proposal of Guth S) which 
would lead to a short period of quasi-exponential expan- 
sion. Amongst cosmologists, the standard wisdom is that 
the classical theory of scalar fields provides an adequate 
description of the relevant physics, and that the new in- 
flationary scenario does not yield an acceptable cosmol- 
ogy. However, a fully quantum-ficld-thcoretic treatment 
even of the simplest model of a phase transition in an ex- 
panding universe is rather difficult, and has rarely been 
attempted. (A series of papers by Boyanovsky, Holman, 
de Vega and others (see Ref. Q and references therein) 
gives a detailed treatment of a model which can essen- 
tially be solved exactly, and I will later discuss this work 
in the context of the programme described here.) 

The work on which I want to give a progress report 
has as its goal the solution of the field equations of semi- 
classical general relativity, namely 



(1.1) 



where Q^ u is the Einstein curvature tensor (which de- 
pends on the metric g and is treated classically) and T M „ 
is the stress tensor of an appropriate quantum field the- 
ory. The simplest version of this problem (which is all 
that I am competent to discuss) takes the spacetime ge- 
ometry to be described by a Robertson- Walker metric, 
with line element ds 2 = a 2 (t) [dt 2 — dx 2 ~\ . Here, t is con- 
formal time, x denotes comoving spatial coordinates and 
a(t) is the scale factor. The field theory I will consider is 
defined in Minkowski spacetime by the Lagrangian 
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C = -d^d^ + -m 2 ^ 2 - -$ 4 , (1.2) 

where $ is a single scalar field, and the bare mass pa- 
rameter TOq is positive, so that spontaneous symme- 



try breaking is expected at zero temperature. In the 
Robertson- Walker spacetime, I suppose this field to be 
conformally coupled to gravity, in which case the scaled 
field 4>(x,t) = a(t)$(x,t) has the action 



d 4 x 



(1.3) 



which is equivalent to a Minkowski-space field theory 
with an effective time-dependent mass m 2 (t) = a 2 (t)mQ. 

The Hamiltonian of this theory is, of course, explicitly 
time dependent. I will take expectation values such as 
that on the right of (1.1) to be 



(0(t)) = Tr [e-* ff (°>0(f)| /Tr [e~' 3off(0) ] , (1.4) 

where 0(t) is a Heisenberg-picture operator whose expec- 
tation value we want to know at time t. I am therefore 
assuming the existence, at some initial time that I call 
t = 0, of an instantaneous state of thermal equilibrium 
at a temperature To = l//3o- Naturally, the initial tem- 
perature is supposed to be high enough for the system to 
start in its symmetric state. Clearly, such an initial state 
is somewhat artificial. I choose it so that I have, at least, 
a well-posed problem that I can attempt to solve. 

I will admit straight away that I have not solved this 
problem. What I will discuss in this talk is, specifically, 
how one may set about estimating the required expec- 
tation values in perturbation theory. For various rea- 
sons, perturbation theory is a rather limited tool for this 
purpose, but I know of no better approximation scheme. 
(The analysis given in [Q avoids the use of perturbation 
theory by restricting itself to the large- A limit of 0(A)- 
symmetric scalar field theory, but does not appear to 
provide a useful basis for a non-perturbative treatment 
of more general models; the functional Schrodinger ap- 
proach described in J5| relies on a Gaussian anzatz for the 
wave functional which seems very difficult to improve on 
in any systematic way and, as far as I know, has been 
made to work only in a universe of 2 spatial dimensions.) 
In setting about this enterprise, I have encountered two 
significant problems. The first is that the propagators 
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of perturbation theory as ordinarily conceived involve 
occupation numbers for quasiparticle modes which do 
not evolve with time as they should when the system 
is driven away from thermal equilibrium. The second is 
that time evolution governed by a symmetrical Hamilto- 
nian (here, H(—(j>) = H{4>)) from a symmetrical initial 
state (in which (<j))=0) cannot produce a nonzero expec- 
tation value for 4>. I will describe how I think these prob- 
lems can be overcome, and how a complete set of Feyn- 
man rules suitable for computing the time-dependent ex- 
pectation values can be obtained. 



II. DISSIPATIVE PERTURBATION THEORY 

My calculations of non-equilibrium expectation values 
are based on the closed-time-path formalism p|-[lO| . More 
specifically, I use the path-integral technique described 
by Semenoff and Weiss fll|| , in which Green's functions 
are obtained from the generating functional 



z (ja) = / [d(j>a\ exp 



iS((/) a ) + i / d xj ■ <f> 



(2.1) 



Here, the single quantum field <p(x,t) is represented by 
three path integration variables <j) a {x, t) (a = 1, ■ • • , 3), 
which live on the Keldysh contour in the complex time 
plane. The action S appearing in Z(j a ) is 



§{<f> a ) = / d 3 x 



dtC( 



dtC((h 



0- 

i / dr Ce{4>z) 
Jo 



(2.2) 



where £(</>) is the orig inal Lagrangian density (in our 
case, that given in (1.2)), while Ce is the Euclidean ver- 
sion associated with the density operator, namely 



4! 



(2.3) 



The time tf is the latest time that we wish to consider, 
and might well be taken to infinity. I shall be particularly 
concerned with the real-time 2-point functions (a, (3 = 
1,2) given by 



G a p(x,x') = 







d 



dj a (x) djp{x' 

Tr[pT(4>(x)$(x'))} 
Tr[p$(x)4>(x>)] 



lnZ(j„) 



Ja=0 

Tr[p$(x')$(x)\ 
Tr[pT(4>(x)4>(x'))] . 

(2.4) 



where T and T denote respectively time-ordered and 
anti-time-ordered products of the quantum field opera- 
tor 4>(x) and p is the initial density operator. I shall 



usually work with the spatial Fourier transforms of these 
functions, which can be written as 

G aP (t, t'; k) = Hp(t, t'; k)6(t - t') + H Q (t' , t; k)6(t' - t) , 

(2.5) 

in terms of a single complex function H(t,t']k), with 
Hi = H and H2 = H*. Other expectation values can, 
of course, be obtained from appropriate derivatives of 
Z(j a ). 

The usual perturbative strategy of taking the quad- 
ratic part of S((p a ) to define the lowest-order approxima- 
tion leads to propagators which can be expressed in the 
same way as the full two-point functions, 

g a p(t, t'; k) = h p (t, t'; k)6(t - t') + h a (t', t; k)9(t' - t) . 

(2.6) 

In this case, the function hk(t,t') is 

h(t,t';k) = {l + n k )f k {t)r k {t')+n k fmh{t') , (2.7) 
where f k (t) is a mode function satisfying 

[d 2 t +k 2 -m 2 (t)\ f k (t) = (2.8) 

and normalized by the Wronskian condition f(t)f*(t) — 
f(t)f*(t) — —i, while n k is the Bose- Einstein occupa- 
tion number n k = (e /3 ° li ' fc — l) associated with the ini- 
tial state. This is clearly unsatisfactory. In the first 
place, since m 2 (t) is positive and increasing with time, 
modes with k 2 < m 2 (t) grow uncontrollably. While such 
behaviour might be expected during the course of the 
phase transition, it cannot be correct in the initial high- 
temperature state. Moreover, the occupation numbers 
are fixed at their initial values, and cannot adequately de- 
scribe the evolving state, except when t is much smaller 
than some characteristic relaxation time. This g a /3 is 
therefore not a good approximation to the full two-point 
function. 

To improve matters, we need mode functions and oc- 
cupation numbers which, to a reasonable approximation, 
describe the elementary quasiparticle excitations of the 
state which actually exists at a given time. To this end, 
I use a more sophisticated means of splitting the action 
into an unperturbed part and an interaction |12-ll4]: 



S(<f>a) = So((f> a ) + Sint(<j)a) 
5*0 W>a) =£ (2) (0a) 



(2.9) 



If f d 3 k 
+ 2 J dtj (2.10) 

Sint=S (4) (0a) 

-\J dtj ^Mk)Ma(i(k)Mk), (2-H) 

where and are respectively the quadratic and 
quartic parts of S. The counterterm J\4 a p is to be 
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chosen, according to an appropriate renormalization 
prescription, to optimize g a p as an approximation to 
G a (3- The new unperturbed action can be expressed as 
So = — i J dbdk(j) a 'D a b4>b-, where I? is a differential op- 
erator, and choosing M. is clearly equivalent to choos- 
ing T>. The real-time propagators g a p are solutions of 
T> ai (t, k)g-yp(t, t'; k) = —iS(t — t'), and it is not hard to 
show that the most general T> which admits a solution of 
the form (EJ]) is 



V, 



a/3 



d 2 + fa - ia k 
-Jkdt 



ia k 



lkd t - 
~d 2 - fa 



- ia k 



(2.12) 



where a k (t), fa{t) and ~f k (t) are undetermined, real 
functions. These functions can be determined self- 
consistently by using the counterterm M. in S- lr ± to can- 
cel some part of the loop corrections to the full two-point 
functions G a p. Details of the prescription used to do this 
are to some extent a matter of choice, and the choices we 
make determine the extent to which properties of G a p 
arc reflected in g a p- For simplicity, I choose P k (f) to have 
the form fa{t) = /x 2 (i) + k 2 , so that fi(t) will be an effec- 
tive quasiparticle mass, but in principle a more accurate 
representation of the quasiparticle dispersion relation is 
possible. The structure of the unperturbed action S ob- 
tained in this way is analogous to that of the effective 
action obtained, for example, in |l5| ] by integrating out 
extra environmental degrees of freedom. 



The propagator is now of the form (2.6), with 



/i(i,f / ;fc) = |exp(-i f t ,dt" lk (t"j) 

x [{N k {t>) + l)/ fc (t)/ fc *(f) + (NW) - l)fmfk(t')\ ■ 

(2.13) 

The mode function f k (t) is a complex solution of 

[d 2 t+ f{t)~\ 1 2 {t) + k 2 ]f k {t) = 0, (2.14) 
while the new function N k (t) satisfies 



d t + 2ift k (t) 



n k (t) 



•7fc(*) 



lk(t)]N k (t) 

= 2ia k (t), (2.15) 



where Q k (t) = l/2f k (t)f^(t). Evidently, the func- 
tion 7fc(£) plays the role of a quasiparticle damping 
rate. When the evolution with time is sufficiently slow, 
and the coupling is sufficiently weak, it is possible to 
identify a time-depende nt o ccupation number n k (t) 



i [N k (t) — 1], and then ( 2.15 ) can be cast approximately 



in the standard form of a Boltzmann equation |13| . The 
dissipative formalism outlined here can be extended to 
complex scalar fields p"6| ] and (with somewhat greater 
difficulty) to spin-i fermions |17|| . 

It will soon be useful to know that, associated with 
T> a/ 3 is an operator 



d a p(t,k) 



which has the property 



d t lk{t) 
-7*(*) ~dt 



(2.16) 



g(t,t";k)d(t",k)g(t",t';k) 

r -ig(t,t';k), t > t" > t' 

= 1 ig(t,t';k), t'>t">t (2.17) 

I , otherwise . 



III. SPONTANEOUSLY UNBROKEN 
SYMMETRY 

As I observed earlier, a state which is initially symmet- 
ric under the operation cj> — > — (j) necessarily evolves into 
another state with the same symmetry if time evolution 
is governed by a Hamiltonian with the same symmetry. 
Nevertheless, as the universe expands and cools, we ex- 
pect to encounter states in which the phenomena conven- 
tionally associated with spontaneously broken symmetry 
(that is, with a non-zero value of ((j))) are realised, and 
I shall refer to such states as having "spontaneously un- 
broken symmetry" . Here, I review briefly the means I 
have found of treating a vacuum state of this kind in 
Minkowski spacetimc |18| . How such a state may be seen 
to emerge from a phase transition is the central question 
that I shall address later on. 

Intuitively, we can imagine a probability density for the 
field 4>(x) at each spatial point x, which has symmetri- 
cally placed peaks both at (f> — er/vA and at cf> = — cr/VA, 
where a j\f\ (with a of order 1) is the expectation value 
normally assigned to in the broken-symmetry vacuum. 
In a large spatial volume, interference between the two 
components of the wave functional is negligible, and the 
doubly-peaked state (with spontaneously unbroken sym- 
metry) should be indistinguishable from the conventional 
broken-symmetry state in which only one of the peaks is 
present. 

More formally, we need a version of perturbation the- 
ory in which the unperturbed Lagrangian describes free- 
particle excitations of the symmetrical vacuum, rather 
than the asymmetrical one. To this end, I take the state 
of spontaneously unbroken symmetry to be characterised 
by the expectation value of </> 2 , which is unconstrained 
by symmetry. Instead of the conventional field variable 
ip(x), defined by 



4>{x) = a + i/j(x) , 
I deal with a field ((x) defined by 

(f> 2 (x) = U 2 + 2U((x), 



(3.1) 



(3.2) 



where U 2 — {4> 2 {x)). Taking m 2 (t) = m§ in (1.3), I find 
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+U (ml ^ ( + l -5\Q) In (1 + 2C/- 1 C) , 

(3.3) 

where the last term comes from the functional Jacobian 
of the transformation and an irrelevant constant has been 
dropped. Since (£) = 0, the linear term must vanish to 
leading order, so we identify U — yM3m^/A(l + 0(A)). 



The expansion of (3.3) in powers of A is 



1 



1 



C=-d^--(2rr4)Q 



(3.4) 



and we see at leading order that particles created by Q 
from the spontaneously unbroken vacuum have the same 
mass V2wo as those created from the broken-symmetry 
vacuum by ip. It is not hard to convince oneself Jl8| that 
these two types of particle arc completely indistinguish- 
able: they have the same mass (located by the pole of 
the propagator) and the same S-matrix elements to all 
orders of perturbation theory. In this formulation, in- 
teractions arise from the kinetic term, and there are an 
infinite number of vertices (though only a finite number 
of these appear at a given order in A). The analysis of 
spontaneously unbroken symmetry outlined here can be 
extended to the Higgs sector of a spontaneously broken 
gauge theory . In accordance with a well-known the- 
orem of S. Elitzur f2(i|l , non-zero expectation values need 
be assigned only to operators which are invariant under 
the gauge and global symmetries of the theory. 



IV. FEYNMAN RULES FOR THE PHASE 
TRANSITION 

As the universe cools, it seems intuitively clear that 
the probability density for 4> starts off with a single peak 
at <p = 0, and subsequently, say at a time to, develops 
the double-peaked form I described earlier. This picture 
is no more that an intuitive guide, because I know of no 
way of actually computing this probability density in per- 
turbation theory. Indeed, the notion that it simply has 
either one or two symmetrically placed peaks may be too 
naive. (It would certainly be too naive to suppose that 
the probability density for values of gives an adequate 
characterization of the total state: for that one would 
need the probability density for functions <j){x) over an 
entire spatial section of the universe, which is much more 
complicated than the one I have in mind.) Specifically, 
the idea I wish to develop is that at times earlier than 
to (whenever this might be), the most probable value of 
</> is zero, and the appropriate way of attempting to use 
perturbation theory is by using the Feynman rules based 
directly on the original Lagrangian ( |1.3] ) in terms of cf> 
itself, while at later times, the most probable values of 



(f> are ±er/ V~\ and the most appropriate perturbation ex- 
pansion is that based on the formulation in terms of £. 
Bringing this idea to fruition is awkward for two reasons. 
One is that the expectation value of <fi 2 is of order 1 in 
the expansion based on 0, while it is of order A -1 in the 
expansion based on £ (or, for that matter, the expansion 
based on ?/>, if there were genuine symmetry breaking). 
Thus, although both expansions are in powers of A, they 
are actually different expansions, corresponding to per- 
turbations about different lowest-order states, and one 
cannot expect to obtain a completely smooth matching 
of these two approximations. My strategy for dealing 
with this is the following. To one-loop order, at least, 
it is possible to identify by eye in the £ formulation a 
quantity a which it seems possible to interpret in the 
way I have described. It should also be possible (though 
I have actually done this only to lowest order) to derive 
an equation of motion for a, roughly of the form 

a(t) = f (a) + A/x(a) + X 2 f 2 (a) + ■■■ , (4.1) 

which is a systematic expansion in powers of A, treating 
a as being of order 1 (though the right hand side is not 
a simple function of <r(t), but rather a functional which 
depends on the values of a at all times between to and 
t). The expectation values we want to calculate can be 
similarly expressed. The initial values of a and a at to 
must, however, be obtained from the values of ((f) 2 ) and 
d((j) 2 }/dt at to as calculated within the expansion based 
on <ft, and are of order yf\. With suc h initial conditions, 
there is certainly a solution to (4.1), but this solution 
cannot be expressed as a power series in A. In any case, 
only a numerical solution is likely to be feasi ble. My 
proposed strategy, then, is to truncate both (hi) and 
the expressions for the initial data at some chosen (in 
practice, no doubt, low) order, and solve the resulting 
problem numerically. 

The second source of difficulty is that the calculations 
using £ for times after to cannot be carried out in isolation 
from those using <f> before to , because the nonequilibrium 
density matrix at to is not directly known. It is there- 
fore necessary to treat the whole problem, making the 
change of variable (3.2) only on the real-time part of the 
closed time path after to- Thus, I need to evaluate a path 
integral of the form 



Z(jaJa) = / [rf0]exp 



i / dt (£(</>)+ j a <f> a ) 



x / [d(] x exp 



i I dt (£(C) + l a <j> a ) 

(4.2) 



This, of course, is not simply the product of two inde- 
pendent path integrals, because of the condition (j> 2 (to) = 
U 2 (to) + 2U(to)((to). Evaluating it is no easy matter, and 
I have not been able to do it as rigorously as I would like, 
but what I find is the following. First, the interactions 
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can be extracted in the usual way, in terms of derivatives 
with respect to the sources: 



Z{ja,l a ) = cxp 



iSt 
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int 

xZ (j a ,l a ). 



exp 



- l Jl 



(4.3) 



The exponential operators of course supply the vertices 
of Feynman diagrams, while the remaining path integral 
Zo(jaja) supplies the propagators. After doing the £ 
integral, I find 



Z = exp 



\ f' dtdt'l(t)g((t,t')l(t') 

Z J to 



to 

exp 



to 

I I [CQ{4>)+ja4>a + (f Z a L a {t) 



(4.4) 



where L{t) is a distribution, concentrated at t = t , such 
that 



to 



dt<j> 2 a {t)L a (t) 



(4.5) 



t=t 



The action for <j> clearly involves a differential operator 
T> a b{l) depending on the source for Q. Apart from some 
details which I will ignore, the result of doing the final 
integration is (in an abbreviated notation which I hope 
is obvious) 

Z = exp [-\lgil - \jg*{l)j - ±Tr In (X>*(i))] . (4.6) 

Expanding g^(t,t'; I) in powers of I, one obtains an infi- 
nite set of propagators. The first is the standard prop- 
agator for (j>, say g^(t,t'). The next is a composite ob- 
ject, consisting of a vertex (involving d a p(t )) located 
at to, from which two <fr propagators and one £ propa- 
gator emerge, and the remainder contain strings of such 
vertices. Similarly, expanding Trln2?(Z) yields a second 
set of composite propagators, consisting of single loops of 
(f> propagators with ( propagators emerging, again from 
vertices located at t . The connection between times 
before t and times after t thus arises from these com- 
posite propagators rather than from vertices. Although 
the propagators are infinite in number, only as many of 
them appear in a diagram of order A™ as can connect the 
vertices which yield this power of A. 

The interactions contained in S? nt naturally consist of 
an infinite set of vertices, just as in the Minkowski-space 
theory. Now, however, the quantity U 2 (t) = {<p 2 {x,t)) 
depends on time. It is still independent of x, owing to 
the translational invariance of the Lagrangian and of the 
initial state. Defining v(t) = A 1 / 2 J7(t), the closed-time- 
path Lagrangian for £ is 



-A- 1 / 2 (v- m 2 (t)v + ^)(( 1 ~C 2 ) 
+ i v -i(v + M 2 (t)v+lv*) (C? - CI) 



(4.7) 



where M(t) is a renormalized mass for the £ excitations. 
The first term constitutes £q, and yields the ( propaga- 
tor. The other two terms that I have shown explicitly 
are counterterms. The one linear in £ can be used to 
enforce the condition (Q = 0, which gives an equation of 
the form v — m 2 (t)v + ^v 3 = loops, to be solved for v(t). 
The term quadratic in £ can be used to cancel an appro- 
priate part of the loop contributions to the £ two-point 
function. Normally, this would constitute a mass renor- 
malization condition. In this formulation of the theory, 
however, the renormalized mass M(t) actually appears 
in the linear term. Nevertheless, these two conditions do 
serve together to give an equation of motion for v and 
a relation between the renormalized mass M(t) and the 
bare mass m(t). 

The equation of motion is somewhat complicated when 
written in terms of v, but simplifies considerably when 
expressed in terms of the quantity 



a(t) = v(t) 



1 - 



2v 2 (t) 



+ 0(A 2 ) 



(4- 



where if (t) denotes the single-propagator bubble 



(4.9) 



In the Minkowski-space theory, Green's functions involv- 
ing the operator tp 2 are rendered ultraviolet finite by a 
suitable combination of additive and multiplicative renor- 
malizations. The same procedure works here (to one-loop 
order, at least) and leads to a renormalized equation of 
motion which is 



a R + M 2 (t)a R - -4 + 0(A 2 ) = . 



(4.10) 



The relation between M(t) and m(t) can also be renor- 
malized, leading to a gap equation 

M 2 {t) = \a%{t) + \ii{t) 

-\a 2 (t)m 2 jl + -A_ [ c - In (a(t)m)} 



+0(A 2 ) , 



(4.11) 



where m is the physical particle mass in the Minkowski- 
space vacuum (which locates the pole of the £ propaga- 
tor) and c = 1 + In 2 - V3n/2. The integral is 
a subtracted version of /£(<), which is ultraviolet finite, 
provided that the large-fc behaviour of the propagators is 
similar to that in the equilibrium theory. In the (/>-based 
perturbation theory before to, one can calculate directly 
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v 2 (t) = I?(t) + 0(\) 



(4.12) 



N k (t) 



2n k (t) + l 



and obtain a renormalized gap equation for the mass fi(t) 
of (j) quasiparticles, which reads 



M 2 (^) = \lf{t) ~ \a 2 {t)m 2 {l + -A- [ c - In (a(t)m)] 



-0(X 2 



(4.13) 



Initial conditions for the solution of (|4.10|) are obtained 
by requiring continuity of v 2 (t) and its time derivative 
at t . The differential equation to be solved for g^(t,t') 
also requires initial conditions. These can be found by 
noting that the connected two-point function for <fi 2 can 
be expressed as 

{^{x,t)4> 2 {x',t')) c = {<t> 2 ( Xl t)ax'X)) c ■ 2\- 1 ' 2 v(t') , 

(4.14) 

and so on, depending on whether its time arguments are 
greater or smaller than to- Initial conditions for g^{t,t') 
are then obtained by requiring this two-point function 
and its first derivatives to be continuous at to- 

In summary, the evaluation of an expectation value in 
the </>-based perturbation theory for early times works as 
follows. We first calculate the Feynman diagrams which 
appear at the desired order in A, obtaining a formal ex- 
pression in terms of the propagators g a (3(t, t';k). In order 
to convert this formal expression into a concrete one, we 
need values for five quantities, namely the mode func- 
tions fk(t), the generalised occupation numbers Nk(t), 
the quasiparticle mass /i(t) and damping rate 7fc(t) and 
the auxiliary functions ak{t). These values are found 
by solving self-consist ently a set of five equations. They 
are: (i) the e quati on ( |2.14 ) for the mode functions; (ii) 
t he eq uation ( 2.15 ) for the Nk(t); (iii) the gap equation 
(4.13) for the quasiparticle mass; (iv) a constraint equa- 
tion for 7/c(t) and (v) a constraint equation for a.k{t)- 
Equations (iv) and (v) are somewhat complicated, so I 
have not given them explicitly (but details can be found 
injOf). They are of the same kind as the gap equation 
( |4.13|) , since all three equations amount to the require- 
ment that u(t), a.k(t) and jk(t) approximate appropriate 
parts of the full self-energy matrix. The equations for 
ctk (t) and 7^ (t) receive contributions from the absorptive 
parts of the self-energy (which are of at least two-loop or- 
der) in place of the one-loop integral if that appears in 
( 4.13 ). Of course, these integrals involve the very quan- 
tities we need to find so, although the solutions are fully 
determined, they cannot be obtained in closed form (ex- 
cept, perhaps, in very special circumstances). In princi- 
ple, however, the equations are susceptible of numerical 
solution. 

Just how difficult this numerical problem is naturally 
depends on what further approximations one is willing to 
make. To give an idea of what is involved, let me consider 
the limit of slow evolution and weak coupling. In that 
case, I can define occupation numbers rifc(i) by 



l- i7fe (t)/2^ 2 +/i 2 W 



(4.15) 



and (|2.15| ) reduces to 



dn 



ak(t) 



Vk 2 + V 2 (t) 



lk{t)n k {t). (4.16) 



In turn, a k {t) and 7^ (t) are giv en by scattering integrals 
which involve nk(t), such that ( 4.16 ) takes the form of a 
Boltzmann equa tion. A t this point , sim ultaneous numer- 
ical solution of ( 2.14 ), ( 4.13 ) and ( 4. 16 ) yields estimates 
of nk(t) and /x(i), from which ak(t) and "fk(t) can be 
reconstructed. The possibility of expressing these two 
functions explicitly in terms of the rik (t) results from the 
fact that they are of order A 2 and can reasonably be set 
to zero within scattering integrals that are already mul- 
tiplied by A 2 . Similarly, 7^(t) might be omitted from 
( [2.14 ) . If one is not satisfied with approximations of this 
sort, then all five equations must be numerically evolved 
together. Two of the equations we have to solve, namely 
( p. 14 ) and ( 2.15| ) are differential equations which need 
initial conditions. These, of course, are obtained from 
the initial equilibrium state and the considerations in- 
volved are discussed in detail in Jl4) . 

Expectation values at later times are obtained from the 
C-based Feynman rules in almost the same way. The dif- 
ferences are: (a) we use the gap equation ( 4.11 ) in place 
of ( |4.13| ); (b) these rules contain a sixth unknown func- 
tion, v(t) (or cr(t)), and, correspon dingl y, there is a sixth 
equation, the equation of motion (4.10), to be solved si- 
multaneously with the other five. Initial conditions for 
this set of equations are obtained from continuity require- 
ments at t = to- 



V. DISCUSSION 

I have explained in outline how I think that pertur- 
bation theory can be used to describe the noncquilib- 
rium dynamics of a symmetry-breaking phase transition 
in an expanding universe. A straightforward application 
of standard perturbative procedures is inadequate for two 
reasons. One is that low-order calculations (which are all 
that one can contemplate tackling) will reflect the evolv- 
ing state of the system only if the propagators are defined 
self-consistently in terms of quasiparticle modes which, 
to a reasonable approximation, represent the elementary 
excitations of the actual state. The other is that sponta- 
neous symmetry breaking as conventionally understood 
cannot come about if the system starts from a symmetri- 
cal state and its evolution is governed by a symmetrical 
Hamiltonian. The Feynman rules I have outlined are de- 
signed to overcome both of these difficulties. 

There are, however, other problems inherent in any 
attempt to treat this problem perturbatively, that I do 
not think can altogether be overcome. One is that the 
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perturbative expansions appropriate to the initial sym- 
metrical state and to the emerging state that I describe 
as having spontaneously unbroken symmetry are really 
two different approximations. This leads to some formal 
difficulties, beyond the ones I have mentioned, that I do 
not propose to discuss here (though I do discuss them in 
a more detailed forthcoming paper) 



no other approach which has this potential. 



A second problem is apparent from equation (4.11), 
where we see a negative contribution to the quasiparti- 
cle mass M 2 (t), whose magnitude increases with time. 
The mode functions from which the pr opag ators are con- 
structed obey equations of the form (2.14) containing a 
squared mass which will be negative for some period of 
time around to, and will grow, potentially to very large 
values. There is clearly the danger that perturbation 
theory will become meaningless if this growth continues 
for too long. The growth cannot continue ind efinit ely, 
because it causes the positive contributions to (4.11) to 
increase in magnitude faster than the negative ones. It is 
therefore not inevitable that perturbation theory need 
become completely unreliable, but concrete numerical 
calculations (which I have not yet undertaken) would be 
needed to assess the extent of the danger. Because of 
this difficulty, Boyanovsky et. al. Q have invested con- 
siderable effort in obtaining numerical solutions for the 
case of the large- N limit of the 0(7V)-symmetric theory. 
The equations they solve (which are exact in the large - 
N li mit) i nclude ones which are rather similar to ( 4.10| ) 
and (4.11), though their meanings are slightly different. 
In particular, a corresponds in that case simply to the 
expectation value of one component of <j), which is identi- 
cally zero for a symmetric initial state. In my language, 
therefore, their calculation is equivalent to working only 
with the <p formulation and this, of course, is possible 
only with an exactly soluble model. The initial condi- 
tions adopted by these authors correspond essentially to 
a vac uum state, with h(t, t'; k) = in place of 

(2.13). While initial conditions appropriate to an initial 
high-temperature state could no doubt also be treated 
by their methods, the corresponding occupation num- 
ber would not be seen to evolve with time, because the 
large- iV theory contains no scattering mechanism for the 
redistribution of energy. For this reason, the final state 
in this calculation is one that oscillates perpetually: the 
thermalization that one would expect in a more realistic 
model does not occur. 

Despite the difficulties I have alluded to, the analy- 
sis I have summarized does lead, at the lowest nontriv- 
ial order, to a well defined set of equations which might 
be solved numerically and should provide a description 
of the quantum dynamics of a phase transition with at 
least the qualitative features (including thermalization 
of the final state) that one would expect. In principle 
(though perhaps with some technical difficulty) the ap- 
proximation scheme can be extended to higher orders, 
and it should also be possible to extend this scheme to 
more realistic models of particle physics. The inherent 
deficiencies of perturbation notwithstanding, I know of 
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